Reversal of the Charge Transfer between Host and Dopant Atoms in Semiconductor 

Nanocrystals 
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We present ab initio density functional calculations that show P (Al) dopant atoms in small 
hydrogen-terminated Si crystals to be negatively (positively) charged. These signs of the dopant 
charges are reversed relative to the same dopants in bulk Si. We predict this novel reversal of the 
dopant charge (and electronic character of the doping) to occur at crystal sizes of order 100 Si atoms. 
We explain it as a result of competition between fundamental principles governing charge transfer in 
bulk semiconductors and molecules and predict it to occur in nanocrystals of most semiconductors. 

PACS numbers: 73.22.-f 
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Introducing appropriate impurity atoms (known as 
"dopants" ) into a semiconductor can dramatically affect 
electrical conduction in the material and is key to the 
operation of modern electronic devices, [lj The dopant 
atoms modify the conductivity of the semiconductor by 
supplying it with additional free electrons or holes that 
can carry an electric current. If an impurity atom hav- 
ing one more electron than an atom of the semiconductor 
host replaces a host atom, in many cases the extra elec- 
tron is very weakly bound to the impurity atom in the 
solid state environment. 2] Thus at room temperature the 
shallow impurity loses (donates) the extra electron to the 
semiconductor and the impurity atom becomes positively 
charged. Conversely an impurity atom with one fewer 
electron than the host accepts an electron from the host, 
becomes negatively charged, and a positively charged free 
hole appears in the semiconductor. This qualitative pic- 
ture of charge transfer between semiconductor host and 
shallow dopant is well established for bulk semiconductor 
materials and is fundamental to our understanding of the 
properties of semiconductor devices. Recent experimen- 
tal and theoretical work HSUBHHUE^has shown 
that it also holds for a variety of doped semiconductor 
nanoparticles. 

Charge transfer also plays an important role in the 
chemistry of molecular systems |l il l, however, the basic 
principles that apply in this case are different: Atoms are 
classified according to their electronegativity which is de- 
fined so that an atom with a larger electronegativity will 
attract (negative) electronic charge from an atom with 
smaller electronegativity. Atoms with nearly filled va- 
lence orbitals have large electronegativities because filled 
orbitals are energetically stable. Conversely atoms with 
nearly empty valence orbitals have low electronegativ- 
ities. Simple semiconductors consisting of atoms from 
group IV of the periodic table have precisely half-filled 
valence orbitals. A group V atom has one more valence 
electron and since its valence orbitals are closer to being 
filled it has a larger electronegativity. Therefore, accord- 
ing to this picture, a dopant atom with one more electron 
than the host semiconductor should attract charge from 
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FIG. 1: On site natural population 19] charge as function of 
the site radial coordinate for Si3sH36. Undoped crystallite in 
top graph. The ellipse in the top graph shows which atoms 
are hydrogen. The central atom (at r — 0) has been replaced 
with aluminum (acceptor in bulk Si) in center graph and with 
phosphorus (donor in bulk Si) in bottom graph. 



the surrounding host atoms and carry a negative charge. 
Similarly a group III dopant atom with one fewer electron 
than the host would be positive. Clearly this reasoning 
based on considerations of molecular chemistry predicts 
charge transfer between dopant and host opposite in sign 
to that found in the solid state semiconductor systems 
discussed above. This raises the intriguing, and until now 
unrecognized, possibility that charge transfer in doped 
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semiconductor nanocrystals with dimensions approach- 
ing the molecular scale might differ fundamentally from 
that in macroscopic semiconductors, with profound im- 
plications for the electronic properties of nanoscale semi- 
conductor devices. Here we show theoretically that this 
is in fact the case and explore the crossover between the 
conventional (macroscopic) and novel nanoscopic/quasi- 
molecular (reversed charge transfer) doping regimes. 

Representative results of our ab initio density func- 
tional calculations ^2 °f the charge distributions in some 
small H-tcrminatcd Si crystals doped with P and Al are 
shown in Fig. ^ Since a P (Al) atom has one more (less) 
valence electron than Si, in bulk Si the shallow dopant P 
(Al) is an electron donor (acceptor) and the impurity site 
is positively (negatively) charged. In Fig. ^ however, the 
reverse is true. Thus our results for these Si nanocrystals 
clearly demonstrate the existence of the quasi-molecular 
regime that we have proposed above, where charge trans- 
fer is governed by electronegativity considerations rather 
than by the standard theory of doping in bulk semicon- 
ductors. 

Ab initio quantum chemistry calculations such as those 
that yielded the results of Fig. Q cannot at present be 
made for much larger crystallites due to practical lim- 
itations of computers. We have therefore developed a 
Poisson-Schrodinger (PS) model for silicon based on a 
nonorthogonal tight-binding (TB) model in order to ex- 
plore the crossover from the quasi-molecular regime to 
crystals large enough that standard semiconductor the- 
ory should be appropriate. We have also examined how 
energy gaps and dopant levels evolve with the size of the 
Si crystallite. 

We have based our TB model on that of Bernstein et 
al.|l3j which reproduces the band structure of silicon very 
well, and gives reasonable values for electron and hole 
masses, see table [I] The on-site potentials in the Bern- 
stein model are however functions of the local density of 
atoms but differences in on-site potentials are explicitly 
given by the PS scheme, so we have used Bernstein's val- 
ues for bulk Si as a starting point for the PS scheme. On- 
site potentials for hydrogen and hopping integrals for Si- 
H have been fitted to reproduce charge distributions ob- 
tained from ab initio density functional calculations. |l2j 
We have used the same overlap and hopping integrals for 
Al-Si and P-Si as for Si-Si. The on-site parameters for 
Al and P have been fitted to yield the correct sign of the 
charge on the impurity site for small crystals and realistic 
values for dopant energy levels for large ones. The on-site 
electron repulsion energies for Al, Si and P are based on 
valence orbital ionization ener gies taken from table D4 of 
Ref. Q; equation (D6) in Ref. Il4l is used for H. We have 
ensured that our model reproduces the on-site energies 
and band structure of the Bernstein model for bulk Si. 



TABLE I: Some properties of the TB model. [l3| experimental 
values are given in parenthesis. 20] 

Position of conduction band minima 87.7% L— X (85%) 

Band gap 1.01 eV (1.12) 

Light hole mass 0.26 m e (0.15) 

Heavy hole mass 0.31 m e (0.54) 

Longitudinal electron mass 0.55 m e (0.92) 

Transverse electron mass 0.15 m e (0.19) 
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FIG. 2: (Color online) The energies of the valence and con- 
duction bands of the crystallites, with and without dopants, 
plotted as functions of crystallite radius. The energies of the 
dopant states are also shown. The inset shows the band gap 
plotted as f(r) = Eg(r)/ E g o — 1 on a log-log-scale for compar- 
ison with papers [l^, [la. A fit,/(r) = 9. Or -10 is also shown. 

The energy on site i includes the electrostatic term 

where e r is the relative dielectric constant due to core 
polarization(core electrons are not included in the TB 
model) in Si, qj is the net charge on site j, is the posi- 
tion of site i, S is the surface of the structure and p(r) is 
the surface polarization charge, due to the core polariza- 
tion. We have chosen the value e r = 6.5 so that the model 
reproduces the correct total dielectric constant of 11.8 for 
infinite Si slabs, taking into account polarization of both 
the valence and core electrons. The electrostatic equation 
is solved self-consistently together with the Schrodinger 
equation Hip = ESip for the nonorthogonal TB model. 
Here S are overlap integrals. The charge on each site is 
calculated using Mulliken population analysis. Q 

We have applied the present model to calculate the 
ground state properties of a number of silicon nanocrys- 
tals ranging in size from from Si2gH36 to Si633H3oo, with 
and without dopants. All of them are approximately 
spherical, have tetrahedral symmetry, and are hydrogen 
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TABLE II: Band gaps for undoped and doped crystallites. 
Dopant levels for the doped nanocrystals. 
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terminated to obtain a clean energy gap for the undoped 
nanocrystals. (ljj 

Figure |3 and Table [H] show how the energies of the 
valence band, conduction band and dopant levels change 
with nanocrystal size. The conduction band energy varies 
little with nanocrystal size and dopant species. The va- 
lence band moves up narrowing the band gap to ~1.7eV 
for Si633H3oo from ~2.7eV for Si29H36. Doping widens 
the band gap somewhat, but this effect is most significant 
for the smallest crystals. The band gap can be fitted to a 
function E g (d)/E g0 - 1 = Ar" 6 , where r = 1.684567V 1 / 3 
is the radius of the crystallite, N is the number of Si 
atoms, E g o is the band gap in the bulk and, A and b are 
fitting parameters. We find that A — 9.0 and b — 1.0; see 
fitted line in inset of Fig. Liu et al.0] and Zunger et 
al.[l6| report b — 1.37 in models without Coulomb inter- 
actions. Effective mass theory (particle in a box) predicts 
an r~ 2 scaling. Our result, however, agrees very well with 
density functional theory calculations 0, 0, Qj| . 

There has been interest in how the system size affects 
the dopant levels and there has been a number of stud- 
ies using different methods, for example effective mass 
theory, [1 |E TB,|j| PRDDO0 and DFT.H3] Our re- 
sults are consistent with this previous work. We define 
the dopant level as the energy difference between the 
partly filled dopant state and its nearest neighbor state. 
The Al dopant levels in the nanocrystals, see Table [HJ 
vary quite smoothly with crystal size (AlSi3iH36 is an 
exception), down to about three times the bulk value 
(57meV) for AlSi632H 30 o- The P dopant levels on the 
other hand vary a lot from cluster to cluster, but also 
reach about three times the bulk value (45 meV) for our 
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FIG. 3: (Color online) Top graph, total Mulliken charge on 
each site as function of radial coordinate for four different P- 
doped clusters. Center graph, Mulliken probability distribu- 
tion for the donor state. Bottom graph, probability distribu- 
tion multiplied by radial coordinate squared to approximate 
probability of finding the electron at a certain radius. 



largest crystallite. We attribute this difference between 
the acceptor and donor states to the fact that the elec- 
tron states have a larger probability on the surface sites, 
making donor states much more sensitive to the surface 
than acceptor states which typically are more localized to 
the interior of the cluster; see the bottom graphs of Figs. 
13 and The strong variations for the donor state levels 
suggest that it would be difficult to engineer the proper- 
ties of a small n-doped cluster without atomic control in 
the manufacturing. It is also relevant in this regard that 
the bands in these small structures are made up from dis- 
crete energy levels and even for our largest crystallites, 
these levels have an energy spacing of 5 — 50 meV. 

The charge on the impurity site (r = 0) for the phos- 
phorus doped crystallites in Fig. |3 exhibits a crossover 
from negative to positive when going from small to 
large nanocrystals. This crossover between the quasi- 
molecular behavior and the bulk semiconductor behav- 
ior occurs between PSi34H36 and PSis6H76. For the alu- 
minum doped crystallites (Fig. we find a crossover 
from a positive to a negative impurity site between 
AlSii22Hioo and AlSii46Hioo- The precise crossover 
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FIG. 4: (Color online) Top graph, total Mulliken charge on 
each site as function of radial coordinate for four different 
Al-doped clusters. Center graph, Mulliken probability dis- 
tribution for the acceptor state. Bottom graph, probability 
distribution multiplied by radial coordinate squared to ap- 
proximate probability of finding the hole at a certain radius. 



points are sensitive to the parameters of the model and 
this result should be regarded as the first (order of magni- 
tude) estimate. The charges on the impurity site in Figs. 
13 and 01 are consistently somewhat smaller in magnitude 
than in Fig. we attribute this difference to the fact 
that Mulliken population analysis (Figs. |3|and0J, tends 
to smear charge between overlapping orbitals on neigh- 
boring atoms more than do natural orbital calculations 

(Fig. m. 

In conclusion: The Poisson-Schrodinger model we have 
developed has allowed us to explore the crossover from a 
novel regime in semiconductor nanocrystals in which the 
molecular view of charge transfer between atoms holds 
true to a regime where macroscopic solid state semicon- 
ductor theory prevails. The crossover is signaled by a 
striking reversal of the sign of the charge transfer be- 
tween the host semiconductor and dopant atom that has 
not been anticipated in previous experimental or theoret- 
ical work. We predict that it should occur at nanocrys- 
tal sizes of order 100 Si atoms. Since very basic prin- 
ciples of solid state semiconductor physics and molecu- 



lar chemistry are the underlying reasons for the charge 
reversal, we predict it to be a general phenomenon oc- 
curring for a wide variety of nanoscopic semiconductors 
and dopants. For Si nanocrystals we also find an energy 
gap widening that scales as r -10 consistent with den- 
sity functional theory calculations. 0, 0, We predict 
the dopant energy levels for Al in Si nanocrystals to vary 
quite smoothly with cluster size while donor levels should 
vary widely from crystallite to crystallite, making it dif- 
ficult to engineer properties of P-doped Si nanocrystals 
without atomic control in manufacturing. 
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